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ABSTRACT 



Aims. To derive the mass profiles of the different luminous and dark components in clusters. 

Methods. The cluster mass profile is determined by using the Jeans equation applied to the projected phase-space distribution of about 3000 
galaxies members of 59 nearby clusters from the ESO Nearby Abell Cluster Survey. The baryonic and subhaloes mass components are 
determined from the galaxies' luminosity-density profiles through scaling relations between luminosities and baryonic and dark halo masses. 
The baryonic mass component associated with the intra-cluster gas is determined using X-ray data from ROSAT. 

Results. The baryon-to-total mass fraction decreases from a value of ^ 0.12 near the center to =^ 0.08 at the distance of =: 0.15 virial radii, then 
it increases again, to reach a value of =^ 0. 14 at the virial radius. Diffuse, cluster-scale dark matter dominates at all radii, but its contribution 
to the total mass content decreases outwards to the virial radius, where the dark matter in subhaloes may contribute up to =^ 23%, and the 
baryons =^ 14% of the total mass. The dark mass and diffuse dark mass profiles are well fit by both cuspy and cored models. The subhalo mass 
distribution is not fit by either model. 

Key words, cosmology: observations - galaxiesxlusters: general - galaxies:kinematics and dynamics - dark matter 
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1. Introduction 



With the increasing accuracy of cosmological numerical 
simulations, the study of the mass profiles of galaxies and 
galaxy clusters has become a powerful way to constrain cos- 
mological models. With numerical simulations, Navarro et al. 
(1996, NFW hereafter) found that dark matter (DM hereafter) 
haloes are characterized by a universal mass-density profile, 
simply summarized by two power-law regimes, an inner one 
with exponent -1, and an outer one with exponent -3. The 
universality of the profile and the existence of a central cusp 
have been confirmed (Moore et al. 1999; Diemand et al. 2004; 
Navarro et al. 2004; but see Ricotti 2003). 

It is obvious that observational knowledge of the mass pro- 
files of the dark and baryonic matter in clusters provides cru- 
cial insights into their formation and evolution (see, e.g., Gao 
et al. 2004; Springel et al. 2001). In particular, one aspect of 
the ACDM theory that can be tested is the presence of a cusp 
in the center of the DM halo. It is well known that on galactic 
scales mass profiles of the NFW form are unable to account for 
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the rotation curves of low-surface brightness galaxies (de Blok 
& Bosma 2002), normal spirals (de Blok et al. 2003; Gentile et 
al. 2004), and the fundamental plane of ellipticals (Borriello et 
al. 2003). 

On cluster scales, the situation is far more open. Cluster 
mass profiles have been obtained from the analyses of the X- 
ray emitting intra-cluster (IC hereafter) gas, of the projected 
phase-space distribution of cluster galaxies, and of the gravi- 
tational lensing shear pattern of background galaxies. Most re- 
sults indicate consistency with the NFW profile (see, e.g., Allen 
et al. 2000; Athreya et al. 2002; Biviano & Girardi 2003; Jee 
et al. 2005; Katgert et al. 2004, hereafter KBM; Lokas et al. 
2006; Pratt & Arnaud 2005; Rines et al. 2003; Rines & Diaferio 
2006). In some cases a flatter than NFW profile is preferred 
(Broadhurst et al. 2005; Ettori et al. 2002; Kelson et al. 2002; 
Nevalainen et al. 2000), or even required (Arieli & Rephaeli 
2003; Demarco et al. 2003; Sand et al. 2004). While isother- 
mal profiles were rejected by some dynamical analyses (e.g. 
Rines et al. 2003; Broadhurst et al. 2005), cored profiles gener- 
ally were not excluded (see however Dahle et al. 2003), as far 
as the core of the matter distribution is small (Arieli & Rephaeli 
2003; Biviano & Girardi 2003; KBM). 

In order to allow a proper comparison of the results of sim- 
ulations with observations, it is important to determine the to- 
tal mass distribution and its major components, baryons and 
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the subhaloes. In particular, the baryonic contribution has been 
shown to be relevant in clusters, both near the centre because of 
the substantial contribution from the cD (e.g. Sand et al. 2004), 
and in the outer regions, because of the increasing mass frac- 
tion of the IC gas (e.g. Lokas & Mamon 2003). Recently, KBM 
have derived a synthetic mass profile of rich galaxy clusters, 
using the ESO Nearby Abell Cluster Survey (ENACS) data-set 
(Katgert et al. 1998). ICBM found that the total cluster mass 
profile is well fitted both by a NEW profile, and by a Burkert 
(1995) profile. KBM also found that the total mass profile is 
very well traced by the luminosity profile of the early-type 
galaxies, i.e. the mass-to-light ratio is almost flat, when only 
early-type galaxies are selected, and the brightest members are 
excluded. The aim of this paper is to derive the mass distri- 
bution of the DM, of the baryonic matter and of the cluster 
sub-haloes, separately. We use the same data-set of KBM anal- 
ysed in a slightly different way (see § I2.1> . We use the X-ray 
data from Reiprich & Bohringer (2002) to determine the ra- 
dial profiles of the IC gas, and the luminosity-density profiles 
of cluster galaxies to determine the mass distributions of the 
galaxy baryons and of the DM subhaloes. 

The structure of this paper is as follows: in §|2we determine 
the mass profiles of the different cluster components; in §|3]we 
fit models to the observed mass profiles; in §|3we summarize 
our results and draw our conclusions. Throughout this paper 
we adopt Hq = 70 km s"' Mpc"', Q.„, = 0.3 and = 0.7. 



2. Mass profiles 

2.1. Total mass 

We use the data-set of 59 clusters used by KBM and described 
in Biviano et al. (2002). These clusters are combined into a sin- 
gle 'ensemble' cluster, in order to improve the rather poor num- 
ber statistics of individual clusters. This gives a total sample of 
~ 2900 member galaxies with positions and redshifts. Note that 
Sanchis et al. (2004) have shown that a stacked sample of sev- 
eral galaxy clusters can be used to determine a reliable average 
mass profile of individual clusters. 

The stacking is done in the space of normalized cluster- 
centric distances, R/rioo, and normalized velocities with re- 
spect to the cluster mean velocity, (v - \')lcry, where rioo is the 
radius of the sphere around the cluster centre with mean den- 
sity equal to 100 times the critical density, and cr,, is the global 
velocity dispersion. Eor the scaling of clustercentric distances, 
Biviano et al. (2002) used Carlberg et al.'s (1997) proxy for r2oo 
(the radius of the sphere around the cluster centre with mean 
density equal to 200 times the critical density), which is based 
on the assumption of an isothermal mass profile. However, the 
cluster mass profile determined by KBM is not isothermal, 
therefore, for the estimation of both rioo and r2oo (needed for 
the determination of the IC gas mass profile, see § l2.2t . we pre- 
fer to use Popesso et al.'s (2005) relation (see eq. 1 in that pa- 
per) which makes explicit use of the shape of the mass density 
profile determined by KBM. The average values of cr,, and rioo 
for our cluster sample are 699 km s"' and 2.25 Mpc, respec- 
tively, and the average virial mass is M(< rioo) = 6.5 10'"* Mq. 
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Fig. 1. The circular velocity profiles, Vc = {GM,o,/rf''^ , of the 
total mass, obtained assuming isotropy (Jl = 0, solid line), 
mild radial anisotropy (yi = 0.1, lower dash-dotted line), and 
mild tangential anisotropy (J?l - -0.6, upper dash-dotted line). 
The dotted fines dehmit the + 1 -cr intervals accounting for both 
the random errors (as obtained from 64 bootstrap resamplings 
of the data in the isotropic case), and systematic errors (arising 
from the uncertainty on the value of the velocity anisotropy y[). 



To determine the mass profile of the ensemble cluster 
(M,o,(r) in the following), we first apply the isotropic Jeans 
equation on the early-type cluster members, as described in 
KBM. The errors on M,o,(r) are obtained from 64 bootstrap 
resamplings of the data (see KBM for more details). Near- 
isotropy for the early-type cluster members was inferred by 
KBM from the analysis of the whole velocity distribution of 
these galaxies. Support for the isotropic assumption also comes 
from the analysis of Biviano & Katgert (2004) of the or- 
bits of different cluster galaxy populations. Here we go be- 
yond KBM's isotropic assumption, and we also consider so- 
lutions with constant orbital anisotropy J?I 0, where = 
1 - <v^>/<v^>, and <Vj^>, <v,^> are the mean squared com- 
ponents of the radial and tangential velocity, respectively (see, 
e.g., Binney & Tremaine 1987). Since KBM constrained the 
velocity anisotropy of the early-type galaxy population to lie in 
the range -0.6 ^ y[ 5 0. 1 , we determine M,o,{r) for the two ex- 
treme cases y[ = -0.6 and ^ = 0.1, using the anisotropic Abel 
inversion (Mamon & Boue 2006) and Jeans equations. The re- 
sults are shown in Eig.^ where we display the circular velocity 
profiles, Vc = {GM,otl rf*'^ ■ While the isotropic solution for the 
total mass profile is not significantly different from that derived 
by KBM, orbital anisotropy has some effect on the uncertainty 
of the total mass profile. 

2.2. Baryonic mass: IC gas 

To determine the mass profile of the IC gas (M^^jsir) in the fol- 
lowing), we use the IC gas density profiles. Unfortunately these 
are available only for a subset of our 59 cluster sample. This 
subset is characterized by a larger average velocity dispersion 
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Fig. 2. The contribution of the IC gas mass to the circular ve- 
locity profile, Vc = (GMgaJr)^'^ (solid line), obtained from the 
sample of Reiprich & Bohringer (2002). The dash-dotted line 
shows the effect of using Neumann's (2005) steeper IC gas den- 
sity profile in the radial range > rsoo- The dashed lines delimit 
the +l-cr intervals accounting for the uncertainties in the best- 
fit parameters of the IC gas density profile, as well as for the 
uncertainty in the value iMga JM,oi){r2oo). 

than the whole sample. Since the shape of the IC gas density 
profile depends on the cluster X-ray temperature (Mohr et al. 
1999), and hence also on the cluster velocity dispersion, we 
cannot use the subset of clusters with available X-ray data as 
representative of the whole sample. 

We proceed instead as follows. We consider Reiprich & 
Bohringer's (2002) sample of clusters for which the best-fit pa- 
rameters of the y6-profiles (Cavahere & Fusco-Femiano, 1978) 
are available, 

PIC ^ poll +(r/ref]-'^"^^^ (1) 

We then consider the distribution of velocity dispersions of 
our 59 clusters, and convert it into a distribution of pseudo X- 
ray temperatures, using the empirical relation of Girardi et al. 
(1996). We generate 500 bootstrap sets of this pseudo-r^ dis- 
tribution. For each of these bootstrap sets, we then extract 59 
clusters from the sample of Reiprich & Bohringer (2002), cho- 
sen to have a r^-distribution as close as possible to the pseudo- 
Tx distribution of the bootstrap set. We finally compute the av- 
erage values of and /? for the 59 x 500 extracted clusters 
(many of the clusters are of course extracted more than once, 
as expected for a bootstrap procedure), < /3 >- 0.625 + 0.007, 
and < Tc >= (0.0636 + 0.0004) r2oo- We use these values to 
define the average gas density profile, which is meant to be 
representative of our cluster sample. Integration of this average 
gas density profile provides the IC gas mass profile, apart from 
a constant, that we fix to the average gas-to-total mass fraction 
at r2oo as determined by Ettori (2003) for a sample of nearby 
clusters, 0.11+J}j;^. 

There is an additional uncertainty that is related to the ex- 
trapolation of the observed IC gas profiles to r2oo- Most of the 



IC gas profiles on which Mgas(r) is based are determined from 
data at radii smaller than ~ rsoo (the radius of the sphere around 
the cluster centre with mean density equal to 500 times the crit- 
ical density) corresponding to ~ 1 Mpc in our cluster sample. In 
order to deal with the systematic uncertainty related to this ex- 
trapolation, we take into consideration the analysis of Neumann 
(2005). Neumann has analysed a sample of 14 nearby clusters, 
for which she has been able to trace the X-ray emission beyond 
r2oo- She concluded that the IC gas density profile steepens with 
increasing radius. Out to 1.2 r2oo Neumann found a best-fit /3- 
profile' with rc = 0.19 r2oo and /3 - 0.8. The fitted profile is 
therefore significantly steeper than the mean profile we have 
found by using the sample of Reiprich & Bohringer (2002). 

We then adopt Neumann's best-fit to estimate the uncer- 
tainty involved in the extrapolation of our average IC gas den- 
sity profile to radii > rsoo- We find that Mgas is not strongly 
modified beyond rjoo (and out to 1 .5 r2oo), despite the substan- 
tial change in the IC density profile, since the normalisation 
of the gas mass profile at r2oo is fixed, iMga JM,oi){r2cn)) - 
0.11 Hence, the uncertainty on the resulting Mgas is mostly 
driven by the uncertainty in the normalization of the gas-to- 
total mass fraction at r2oo- In Fig.|2lwe show the resulting con- 
tribution of the IC gas mass component to the circular velocity 
profile, with its confidence interval, accounting for the uncer- 
tainties in the best-fit parameters of the IC gas density profile, 
and in the normalisation value (Mgas/M,o,)(r2oo). 

2.3. Baryonic mass: galaxies 

In order to determine the baryonic mass profile in the galac- 
tic component {M^^'"^{r) hereafter) we consider separately the 
two classes of early- and late-types (we use the data-set of 
Thomas & Katgert 2006 for the morphological/spectral clas- 
sification of ENACS galaxies). The baryonic mass profiles of 
these two galaxy classes are obtained following the prescrip- 
tions of KBM for the derivation of their number density pro- 
files (see Appendix B in KBM; see also Biviano & Katgert 
2004). KBM's methodology accounts for both the incomplete 
azimuthal coverage of the ENACS observations (Katgert et al. 
1996, 1998), and the fact that different clusters are sampled out 
to different fractions of their virial radii. 

In KBM the number density and the luminosity density pro- 
files were derived; here we proceed further by assuming (for 
each galaxy class) a relation converting a given galaxy lumi- 
nosity to its baryonic component. For the early class we take 
Borriello et al.'s (2003; see their eq.l5) relation, and convert 
their magnitudes to the R-hwd magnitudes of ENACS galax- 
ies (Katgert et al. 1996), using the relations of Fukugita et al. 
(1995). For late-type galaxies the relationship between lumi- 
nosities and baryonic masses is taken from Salucci & Persic 
(1999; see their eqs. 3-5). 

Since the ENACS sample is not complete, we need to cor- 
rect the density profiles for incompleteness. We estimate that 
the ENACS sample is roughly 75% complete down to an appar- 

' Among the subsamples of Neumann (2005) we have chosen her 
no.l, based on the similarity between the Tx distribution of this sub- 
sample and of our clusters. 
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Fig. 3. The contributions to the circular velocity profile of the 
baryonic mass in all cluster galaxies (solid line), in early-type 
cluster galaxies (dash-dotted line), and in late-type galaxies 
(dashed line). The dotted line shows the profile derived by 
Lokas & Mamon (2003) for the Coma cluster, after appropriate 
scaling (see text). 



ent magnitude R - 16.5 (see Fig. 4 in Katgert et al. 1998), and 
then the completeness rapidly drops for fainter magnitudes. At 
the median redshift of our clusters (z - 0.064), R - 16.5 corre- 
sponds to an absolute magnitude Ms = -20.8. For simplicity, 
we then assume 0.75 completeness down to Mr = -20.8, and 
zero at fainter magnitudes. 

In order to derive the luminosity of galaxies fainter than 
M/f = -20.8, we use the 7?-band luminosity function of 
Lugger (1986), that is a Schechter (1976) luminosity function 
with M^ = -21.9 and a - -1.24. We integrate it between 
Mr - -20.8 and the magnitude corresponding to O.OIL^, 
where is the luminosity corresponding to M^. We further 
assume that most faint cluster galaxies with M^ > -20.8 are 
dwarf spheroidals, i.e. early-type galaxies. We find that galax- 
ies fainter than Mr = -20.8 contribute 25% of the galactic 
baryonic mass in a cluster We then correct the observed bary- 
onic profiles for the faint galaxies contribution, and for the ad- 
ditional factor 1 .33 = 1 /0.75, that accounts for the overall spec- 
troscopic incompleteness. 

The resulting contribution of the baryonic mass in cluster 
galaxies to the circular velocity profile is shown in Fig. |3] In 
the same figure we also show the separate contributions of the 
baryonic mass in early-type and late-type galaxies. The profile 
of the baryonic mass in early-type galaxies is similar to that 
of the total mass, but more centrally concentrated. The profile 
of late-type galaxies is instead less centrally concentrated than 
the total mass profile. Early-type galaxies dominate the galactic 
baryonic budget within the virial radius, but the contribution 
of late-type galaxies to this budget increases with radius, from 
only 3% near the center to ~ 40% within the virial radius. 

In order to put our galaxy baryonic mass profile in context, 
we compare it with that derived by Lokas & Mamon (2003) for 
the Coma galaxy cluster. First we compute the ratios between 



the Coma virial radius and total mass and the corresponding 
average values for our 59 clusters (these ratios are 1.3 and 1.9, 
respectively). We then rescale Lokas & Mamon's profile using 
these scaling factors, and assuming that the cluster baryonic 
mass scales proportionally to the cluster total mass. The result 
is shown as a dashed line in Fig.|3l Clearly, Lokas & Mamon's 
profile is less concentrated than ours, but it does not differ by 
more than +20% over the whole radial range considered here. 
We consider this a remarkable agreement, given the fact that 
Coma is not a typical galaxy cluster. While 80% of our clus- 
ters are dominated by a single brightest cluster galaxy (BCG 
hereafter), two BCGs are present in the inner region of Coma, 
which is in fact substructured (see, e.g., Adami et al. 2005 and 
references therein). The presence of substructures in the Coma 
cluster core probably reduces the concentration of its baryonic 
mass profile. 

2.4. Dark mass: subhaloes 

The distinction between DM in subhaloes and diffuse DM is 
to some extent a matter of definition. Here we consider as sub- 
haloes all visible cluster galaxies (including those we miss be- 
cause of the completeness limit of the ENACS survey), with 
the exception of the BCG. Since the BCG is centrally located, 
its diffuse stellar and DM haloes are generally considered part 
of the diffuse IC material (Lin & Mohr 2004; Murante et al. 
2004). 

We compute the subhalo mass profile (Msub(r) in the fol- 
lowing) from the average luminosity density profiles of the 
cluster member galaxies, by adopting a Hubble-type dependent 
scaling relation between a galaxy luminosity and its halo mass 
(see Shankar et al. 2006). As in § 12.31 we apply the needed 
corrections for incompleteness. Galaxies fainter than the spec- 
troscopy limit of the ENACS survey (Mr ^ -20.8) contribute 
roughly 30% of the subhaloes mass. An additional factor of 
1.33 must be included to account for the average incomplete- 
ness of the ENACS spectroscopic sample (see § l2.3t . 

The resulting subhalo mass profile needs however to be 
modified to take into account the fact that, unlike in the field, 
galaxies in clusters are so densely packed that their haloes 
would overlap. This overlap does not in fact occur because their 
haloes are tidally stripped as they pass through dense regions. 
Numerical simulations predict a flattening of the density pro- 
files of subhaloes near the cluster centre, and as much as a 50% 
loss of the total mass in subhaloes, due to the stripping process 
(see, e.g., Gao et al. 2004). Gravitational lensing observations 
confirm that the haloes of galaxies near the cluster centres are 
indeed less massive than those of field galaxies (Natarajan et 
al. 2002; Gavazzi et al. 2004). 

In order to take into account the effects of tidal stripping 
near the cluster centre, and given the current, rather loose, ob- 
servational constraints on this topic, we try the following sim- 
plified approach. We modify our estimate of the subhalo mass 
as follows. First, we assume that subhalo stripping only oc- 
curs within a critical radius Rbcg^ and that stripping is more 
effective as the clustercentric radius decreases. In practice, we 
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Fig. 4. The contribution of the mass in subhaloes to the circular 
velocity profile, Vc = (GM sub I rf'^ (solid line) and its 1-cr con- 
fidence interval (dashed lines) accounting for both the random 
and the systematic errors (the latter being related to the un- 
certainty in the size of the central stripping region). The dash- 
dotted line corresponds to the scenario in which subhaloes lose 
50% of their total mass. 



assume Msub ~ at the cluster centre, and smoothly interpolate 
to the unstripped Msub at the radius Rbcc- 

We consider it reasonable to identify the value of Rbcg with 
the radius of the halo of the BCG. BCGs can be followed pho- 
tometrically out to ~ 0.4 Mpc (see, e.g., Gonzalez et al. 2005). 
The haloes of BCGs have been modelled with a Burkert pro- 
file with scale radius Tj - (2.6 + 0.2)/?f (Borriello et al. 2003), 
where Re - 18+2 kpc (Nelson et al. 2002), implying that 1/2 
of the total BCG mass is contained in a sphere of ~ 0.5 Mpc ra- 
dius. Moreover, numerical simulations predict that the density 
profile of subhaloes deviates from the overall density profile of 
DM inward of ~ 0.4 raoo - 0.7 Mpc (see Fig. 1 in Gao et al. 
2004). Therefore we take Rbcc ~ 0.5 + 0.2 Mpc. The uncer- 
tainty in the value of Rbcc increases the error of M^ub in the 
central region. 

The resulting contribution of M^ub to the circular velocity 
profile is shown in Fig.@] We also display the circular velocity 
profile corresponding to the rather extreme scenario in which 
50% of the total subhaloes mass is lost at any radii {strong 
stripping scenario hereafter.) 



2.5. Dark mass: diffuse 

The diffuse DM profile is obtained from the total mass profile 
by subtracting from it the baryonic mass profile, Mdark(r) = 
M,o,(r) - M^J(r) - Mj„,(r). 

Baryons are a minor though not irrelevant component of the 
total cluster mass. On the other hand, a non-negligible fraction 
of the total cluster mass is in subhaloes. When both the bary- 
onic and subhalo masses are subtracted from the total mass, 
we obtain what we call the diffuse DM profile M'^^Ur) = 



2.6. Relative fractions 

The mass fractions of the different cluster components relative 
to the total cluster mass are shown in Fig.|5]as cumulative mass 
fractions within the radii on the x-axis. For clarity, confidence 
intervals are not shown. We only show the error-bars at r2oo- 
Error intervals at other characteristic radii are listed in Tabled 
(see also Figs.^Q. The largest uncertainty is on M,o,, and this 
affects M^^^.^ (see also Fig.|S|l. 

Among the baryonic components, the IC gas is clearly 
dominant at all radii. However, very near the cluster center, 
galaxies contribute almost as much baryonic mass as the gas. 
This is due to the presence of the cD and/or very bright galax- 
ies near the cluster center (the phenomenon also known as 'lu- 
minosity segregation', see, e.g., Biviano et al. 1992, 2002). 
Baryons in galaxies and the IC gas are more, respectively less, 
centrally concentrated than the total mass; hence, the ratio of 
the total baryonic mass to the total cluster mass has a minimum 
at ~ 0.15 rioo - 0.3-0.4 Mpc. The ratio between the baryonic 
and total mass profiles is constant to within +30% out to the 
cluster virial radius. 

The mass contribution from the subhaloes is small near the 
center, where stripping may occur (see § \2.A\ . and then in- 
creases out to 1 . 1 Mpc, when it reaches an approximately con- 
stant fraction of ~ 0.2-0.25. 

Diffuse DM is the dominant mass component at all radii. Its 
contribution to the total mass almost monotonically decreases 
with radius, except perhaps at the center, where the galactic 
baryons contribute significantly. 

The relative contributions of diffuse and subhaloes DM 
change if one considers the strong stripping scenario. In this 
scenario, the contribution of subhaloes to the total mass budget 
is similar to that of all baryons (see the values listed in brackets 
in TableQ and diffuse DM contributes 3/4 of the total mass. 

3. Model fits 

We fit the observed DM profiles with two models, the NFW and 
Burkert profiles. Both models are characterized by the same 
asymptotic slope at large radii, p{r) cc r^^, but the former is 
characterized by an inner cusp, the latter by a core. By compar- 
ing the results of the two model fits we can address the debated 
issue of the reality of the inner cusp that numerical simulations 
predict to exist in DM haloes. 

The NFW mass density profile model can be written as; 

Po 



Pnfw 



(r/r,)(l + r/r,)2 



and the Burkert profile as: 



PBurkert 



PO 



(l+r/ro)[l+(r/ro)2] 



(2) 



(3) 



MUr) - M'l"l(r) - MgUr) - MM- 



They are both characterized by a scale radius, r s or ro, which, in 
the case of the NFW profile is usually referred to as the inverse 
of the concentration parameter, c = rioo/^s, and, in the case 
of the Burkert profile, corresponds to the radius at which the 
central density drops by a factor 4. 

In order to determine the best fit parameters of the two pro- 
files, we fit the observed circular velocity profiles, Vc- These 
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eiTor bar on {Mg„s/M,o,){r). 



Mass 

component 


0.1 rioo 




fraction of total mass within 

''500 r 200 


'■lOO 




0.07 ± 0.03 




0.09 ± 0.03 


0.11 ±0.03 


0.12 ±0.03 


M'""! 

gal 


0.04 ±0.01 




0.019 ± 0.003 


0.020 ± 0.002 


0.021 ±0.002 




0.04(0.02) + 


0.04 


0.22(0.11) ±0.04 


0.23(0.12) ±0.03 


0.23(0.12) ±0.03 


dork 


0.85(0.86) ± 


0.41 


0.67(0.78) ± 0.20 


0.64(0.75) ±0.14 


0.63(0.74) ±0.13 



Table 1. Relative contributions of the different cluster mass components at four characteristic radii. Values in brackets are for the 
strong stripping scenario. 



are displayed in Fig.|6lfor the different mass components. In 
order to determine the quality of the fit, we make use of the 
standard analysis, applied on the mass density profiles, p(r), 
instead of on the circular velocity profiles, Vdr), since the error 
bars of p{r) at different radii are independent of one another. 



We start by considering the circular velocity profile corre- 
sponding to Mdarkif)- The best-fit NFW and Burkert profiles 
are displayed in Fig. top-left panel. The results are given 
in Table |2] where the characteristic radii are all given in units 
of rioo (- 2.25 Mpc). Both the NFW and Burkert models 
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Fig. 6. The contributions of the different cluster mass components to the circular velocity profile, Vc = (GM/r) . Upper (thick) 



soUd line: M,ou lower (thin) soUd line: M^^-^^ ; dotted line: M^™'; dashed Une: Mjo.,; dot-dashed line: M^™^+Mgas^, triple dot-dashed 
line: Msub- For clarity, only representative 1-cr error bars (at the raoo radius) are displayed. Note that the error bar on the galaxies' 
baryonic profile is very small and not visible in this plot. The error bar on the total baryonic profile is not shown, since it is almost 
identical to the error on the IC gas baryonic profile. 



provide acceptable fits to the circular velocity profile corre- 
sponding to Mdarkif)- The best-fit NFW profile has a a con- 
centration parameter c = 7 + 1, which is in line with the pre- 
dictions of cosmological simulations (c = 8.3 for clusters of 
M(< rioo) = 6.5 10'"* Mq, see Dolag et al. 2004 2). 

To fit the circular velocity profile corresponding to the dif- 
fuse DM component, both the NFW and the Burkert mass 
distributions become more concentrated (see Fig. top-right 
panel, and Table O. This is because the subtracted subhalo 
mass distribution is less concentrated than the total DM distri- 
bution. The best-fit ro value of the Burkert profile is small, but 
still significantly different from zero. The best-fit NFW pro- 
file has a concentration c - lO^j- The value of concentration 
does not change significantly, if we consider the strong strip- 
ping scenario (c = 8 + 2, see Fig. bottom-left panel, and 
values in brackets in Table|2jl. 



^ Dolag et al.'s (2004) definition of c is diff'erent from that used here; 
we converted their values to make them consistent with the definition 
of c used in the present paper. 



Mass profile 


Model 


r, or ro 
(rioo units) 


reduced 




NFW 


0.15 ±0.03 


0.3 




Burkert 


0.09 ± 0.01 


04 


ft 


NFW 


0.10(0.13) ±0.03 


0.3 (0.2) 




Burkert 


0.07(0.08) ± 0.02 


0.3 (0.3) 




NFW 


0.33 ± 0.05 


2.6 




Burkert 


0.19 ±0.02 


2.4 



Table 2. Results of model fits to the observed mass profiles. 
Values in brackets are for the case that 50% of the total mass in 
subhaloes is tidally stripped. 



Neither the NFW nor the Burkert model provides accept- 
able fits to the circular velocity profile corresponding to Mguhir) 
(see Fig.0 bottom-right panel, and Table O. The rather large 
best-fit values of and ro (for the NFW and Burkert model. 
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Fig. 7. Top-left panel: The circular velocity profile, Vc = {GMdarkI rf'^ , corresponding to the cluster DM profile (filled squares 
with 1-cr error bars), and the best-fitting Burkert (solid line) and NFW (dashed hne) profiles. Top-right panel: same as top-left 
panel, for the circular velocity profile corresponding to M^'J/^. Bottom-left panel: same as top-left panel, for the circular velocity 
profile corresponding to the cluster diffuse DM profile, assuming tidal stripping of 50% of the subhaloes mass. Bottom-right 
panel: same as top-left panel, for the circular velocity profile corresponding to Msub- 



respectively) reflect the fact that subhaloes tend to avoid the 
central cluster region (this was already quite clear from Fig.|5}. 

4. Summary and conclusions 

We have obtained the mass profiles of the dififerent cluster com- 
ponents, namely the baryons (in galaxies and the IC gas), the 
subhaloes (galactic halo DM), and the diffuse DM, out to the 
cluster virial radius, rioo - 2.25 Mpc. 

We have determined a total budget of the different mass 
components at different radii. The diffuse, cluster-scale, DM 
is the dominant cluster mass component at all radii. The total 
baryonic mass fraction, resulting from the summed contribu- 
tion of the galactic and IC baryons, is 14% of the total clus- 
ter mass within the virial radius. The baryonic mass fraction 
first decreases, and then increases again with radius, changing 
by +30% within the cluster virial radius. The galaxy baryonic 
component is always a small amount of the total mass, ex- 
cept near the cluster center, where galaxies contribute almost 
as much baryonic mass as the IC gas (because of the centrally 
located cD). Considering clusters as a cosmic laboratory, the 



low relative fraction of baryons in galaxies indicates how the 
star-formation in the Universe has been inefficient, since only 
^ 14% of the baryonic mass content in clusters has been trans- 
formed in long lived stars. 

The baryonic mass profile of early-type galaxies has a 
shape similar to that of the total mass profile (in agreement 
with e.g. van der Marel et al. 2000, KBM, Biviano & Girardi 
2003, Lokas & Mamon 2003). On the other hand, the bary- 
onic mass profile of late-type galaxies is less concentrated than 
the total mass profile. Most galaxy baryons in clusters are con- 
tributed by early-type galaxies. The IC gas-to-total mass frac- 
tion increases with radius as r*'"* beyond ~ 0.2 rioo, in agree- 
ment with previous findings (e.g. Allen et al. 2000, Lokas & 
Mamon 2003). 

We estimate the subhalo mass fraction, 12-23%, which is 
in agreement with the value predicted by the results of cosmo- 
logical numerical simulations (see, e.g. Takahashi et al. 2002; 
van den Bosh et al. 2005; Gill et al. 2004), and estimated ob- 
servationally by Natarajan et al. (2002). 

We obtain good fits to the DM and diffuse DM profiles with 
the NFW and Burkert models, while the subhalo mass profile 
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cannot be fitted by these models. For the DM profile we obtain 

MiarklilO'^ M^) = 4.7 [ln(l + r/0.33) - 

r/0.33(l +r/0.33)"'] (4) 

with r in Mpc, for the NFW fit, and 

MdarkKlO'^ M^) = 1.6{ln(l +r/0.21)- 
arctan(r/0.21) + 0.5 ln[l + (r/0.21)2]) (5) 

for the Burkert fit, for an average cluster of mass M,o,(< rioo) = 
6.5 10"* Mo, and rioo = 2.25 Mpc. For the diffuse DM profile 
we obtain 

^li(/(10"^o) = 2.8 [ln(l + r/0.24) - 

r/0.24(l +r/0.24r'] (6) 

for the NFW fit, and 

m2(/(10'^Mo) = 1.3 {ln(l + r/0.19) - 

arctan(r/0.19) + 0.5 ln[l + (r/0.19)^]) (7) 

for the Burkert fit. In the strong stripping scenario (see § I2.4l i 
we obtain instead 

M^^f/ClO'^Mo) = 3.6 [ln(l + r/0.28) - 

r/0.28(l + r/0.28)-^] (8) 
for the NFW fit, and 

M^^f/dO'^Mo) = 1.3 {ln(l + r/0.19) - 

arctan(r/0.19) + 0.5 ln[l + (rlQ.l9f]} (9) 

for the Burkert fit. 

Note that the best-fit NFW concentration values (c = 7 + 1 
and c - 10^2' respectively, for the DM and DM diffuse profiles) 
are similar to the predictions of ACDM cosmological simula- 
tions (e.g. Dolag et al. 2004). 

Our results are for clusters on average, and it is possible that 
individual clusters are characterized by different types of mass 
profiles (see, e.g., Ettori et al. 2002). Our results provide new 
constraints relative to most recent analyses (e.g. van der Marel 
et al. 2000; KBM) in that we derive the dark, rather than the 
total mass distribution. Lokas & Mamon (2003) did consider 
the dark matter distribution, but only of one cluster (Coma), 
while we have examined a sample of 59 clusters. Moreover, 
they did not subtract the contribution of subhaloes from the 
DM profile. 

The major uncertainty in our analysis is due to the still 
rather limited number of available cluster galaxies with red- 
shifts. Other sources of uncertainties (such as systematics in the 
knowledge of the IC gas density profile at large radii, or of the 
subhalo stripping efficiency), are smaller With a larger data-set 
it will be possible to improve the constraints on the cluster mass 
profile, due to a better definition of the projected phase-space 
distribution of cluster members, and their orbital anisotropy. 
Improvements over the current analysis are therefore expected 
from the use of larger samples of clusters extracted from the 
Sloan Digital Sky Survey (e.g. Abazajian et al. 2004), which 
should allow a factor ~ 5 increase in size with respect to the 
ENACS data-set used here (see Goto 2005). 
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